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How effective are governmental incentives to achieve widespread vaccination coverage so as to prevent 
epidemic outbreak? The answer largely depends on the complex interplay among the type of incentive, 
individual behavioral responses, and the intrinsic epidemic dynamics. By incorporating evolutionary games 
into epidemic dynamics, we investigate the effects of two types of incentives strategies: partial-subsidy policy 
in which certain fraction of the cost of vaccination is offset, and free-subsidy policy in which donees are 
randomly selected and vaccinated at no cost. Through mean-field analysis and computations, we find that, 
under the partial-subsidy pohcy, the vaccination coverage depends monotonically on the sensitivity of 
individuals to payoff difference, but the dependence is non-monotonous for the free-subsidy policy. Due to 
the role models of the donees for relatively irrational individuals and the unchanged strategies of the donees 
for rational individuals, the free-subsidy policy can in general lead to higher vaccination coverage. Our 
findings indicate that any disease-control policy should be exercised with extreme care: its success depends 
on the complex interplay among the intrinsic mathematical rules of epidemic spreading, governmental 
policies, and behavioral responses of individuals. 



Vaccination is one of the most effective tools for reducing morbidity and mortality associated with infectious 
diseases. In spite of the fact that voluntary vaccination may be infeasible to prevent pandemics', global 
vaccination programs for measles, pertussis, and polio have reduced the prevalence of these diseases 
dramatically over the last decades^ '. However, vaccination represents a long-standing social dilemma for public 
health administration, because voluntary vaccination often cannot result in sufficiently strong herd immunity for 
disease eradication'' but compulsory vaccination may result in infringement of civil rights. In general, vaccination 
protects not only those who are vaccinated but also their neighbors. As a result, many others in the community 
can also be benefited, including those who have less incentive to be vaccinated. This scenario naturally leads to the 
"free-riding" problem commonly observed in public goods studies^. However, in voluntary vaccination, an 
individual's decision-making with respect to vaccination may depend on many factors such as the perceived risk 
of infection, cost of infection, cost of vaccination, and the reactions of other vaccinated individuals. Generally, it 
has been demonstrated that a voluntary vaccination policy without incentives can hardly be effective to eradicate 
vaccine-preventable diseases'" '^. A possible resolution of the social dilemma with respect to vaccination is for the 
government or health organizations to implement some incentive programs, e.g., subsidy and insurance policies. 
In this regard, many quantitative studies addressed the effects of the incentive programs on disease eradication, 
and some optimal strategies were proposed'""^". 

There has been a great deal of recent interest in understanding the effects of incentives through mathematical 
modeling and computations. One general finding is that external incentives may alter the vaccination decisions of 
the broader public, thereby reducing the effectiveness of such programs or even making them detrimental. For 
example, by constructing an individual-level adaptive decision-making model based on the minority game 
theory, Vardavas et aU* found that severe epidemics cannot be prevented unless vaccination programs offer 
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incentives. While incentive policies such as offering several years of 
free vaccines to individuals who pay for one year of vaccination can 
substantially reduce the frequency of severe epidemics, incentive 
programs that provide free vaccination to families could increase 
the frequency of severe epidemics^'. In Ref 22, Wells et al. showed 
that the effects of the so-called immunization strategies to vaccinate 
superspreader can be weakened or counteracted when other indivi- 
duals make vaccination decisions about influenza based on the self- 
interested principle. In Ref 23, authors assumed that farmers make 
decisions on whether or not to vaccinate their own herd based on 
self-interested rules such as the payoff maximization rule^"'^* and 
showed that voluntary vaccination cannot eradicate disease among 
herd. Two incentive programs were then considered, namely subsid- 
izing vaccination and compensating losses. It was found that the 
former can lead to a marked suppression of the disease, but cannot 
completely eradicate it. On the contrary, the compensating losses 
strategy has little impacts on the reduction of diseases. Recently, 
the effects of two types of subsidy policies, free- and partial-subsidy 
policies, were studied^^ based on the self-interested rule. The finding 
is that, due to the "herd immunity" effect, the free-subsidy policy on 
promoting the vaccination coverage is not as effective as the partial- 
subsidy policy. 

This work is motivated by the intuition that human behavioral 
responses, naturally, would have some significant impact on the 
complex dynamical interplay among epidemic spreading, vaccina- 
tion, and incentives. When facing the outbreak of an epidemic, it is 
often difficult or even impossible for individuals to act with perfect 
information to assess their chances of becoming infected. In this case, 
individuals may adopt new strategies through learning. The imita- 
tion (learning) rule has been fully studied in evolutionary game 
theory^'""^'', and also been incorporated in many vaccination 
decision-making models, with the general finding that the rule can 
yield results different from those based on the self-interested 
principle'''"''^. 

In this paper, we investigate how the imitation rule affects two 
previously studied subsidy policies: (1) the partial-subsidy policy, 
where a certain proportion of subsidy is distributed to all vaccinated 
individuals and (2) the free-subsidy scenario in which, given the 
limited amount of subsidy, a certain number of individuals are ran- 
domly vaccinated without any personal cost. We focus on a flu-like 
vaccination and develop an epidemiological game-theoretic model 
by integrating a classical epidemic-spreading process and the subsidy 
policy into a simple agent-based model. Our analysis and computa- 
tions indicate that the cost of vaccination and strength selection in 
the imitation rule play a vital role in the effectiveness of the subsidy 
policies. More specifically, for small cost of vaccination, the vaccina- 
tion coverages for the two policies increase with the selection 
strength characterizing the rational/cautious behaviors in the 
imitation rule. On the contrary, for large cost of vaccination, increas- 
ing the value of P gives rise to different results for the two subsidy 
policies: for the partial policy, the vaccination coverage increases 
with ji when the proportion of subsidy is small and the opposite case 
occurs when the proportion of subsidy is considerable. For the free- 
subsidy policy, the vaccination coverage is non-monotonically 
dependent on the value of ji. Owing to the role models of donees, 
the vaccination coverage for the free-subsidy case is generally higher 
than that of the partial-subsidy case. These results provide insights 
into the role of subsidy policies associated with vaccination in pre- 
venting large-scale epidemic outbreak. 

Results 

Model description and analysis are described in Methods. The para- 
meter fi in Eq. (1) is the strength of selection (0 < P < «:), which 
characterizes the sensitivity of individuals in response to payoff dif- 
ferences, where a larger value of P means the individuals are more 
rational (less random) with respect to making decisions about vac- 



cination. We first present results with homogeneous small-world 
networks (HSWN)^''''", where we study the effects of different values 
of /i on the two subsidy policies, as shown in Figs. 1 and 2 for c = 0.5 
and c = 0.9 (we have checked the smaller values of c and found that 
the phenomenon persists. For example, the case of c = 0.2 is given in 
Fig. SI.), respectively, where c is the cost of vaccination [Eq. (1)]. We 
observe distinct effects on the vaccination coverage (V) and the epi- 
demic size (R) for different costs of vaccination. 

For small vaccination cost (c = 0.5, Fig. 1), as fl is increased, the 
vaccination coverages (V) for the partial- and free-subsidy policies 
decrease monotonically and then reach a constant, leading to an 
increase in the epidemic size. For small values of P, the blind imita- 
tion behaviors of individuals in combination with small vaccination 
cost give rise to a relatively high level of vaccination coverage, result- 
ing in a quite low level of the epidemic size. When P becomes larger, 
individuals are more sensitive to the payoff difference. In this case, on 
average those taking risk can get higher payoff than those taking 
vaccination since the epidemic size is small. As a result, more indi- 
viduals tend to imitate the strategy of free riders. Eventually, the 
vaccination coverage is reduced, leading to an increase in the epi- 
demic size. If the vaccination coverage keeps decreasing, then the 
epidemic size would increase continuously, generating greater risk 
for the non- vaccinators. Driven by rational thinking, certain indivi- 
duals would choose to be vaccinated, leading to a turn-around in the 
number of vaccinated individuals. The two competing factors, 
namely higher payoff and risk, eventually make the vaccination cov- 
erage to converge to a stable level for large values of p. The stable, 
equilibrium value of the vaccination coverage can be theoretically 
predicted based on the mean-field approximation [Eq. (18) and Eq. 
(20) in Methods]. The insets of Fig. 1 give the differences in the 
vaccination coverage and epidemic size between the partial- and 
free-subsidy policies, Vp — Vp and Rp — Rp, respectively, obtained 
from direct numerical simulation and theory. We observe reasonable 
agreement between numerics and theory. 

Results for relatively large vaccination cost (c = 0.9) are more 
complicated than those for c = 0.5, as shown in Fig. 2. In particular, 
for the partial-subsidy policy, the vaccination coverage increases with 
P when the subsidy is not substantial [e.g., S = 0.1 in Fig. 2(al) and 
5 = 0.4 in Fig. 2(a2)] , where 5 is the fraction of total vaccination cost 
contributed by subsidy. For small values of p, individuals are not 
quite rational, so they are not willing to take vaccination due to the 
high cost, generating a relatively low level of vaccination coverage 
and, consequently, larger epidemic size. For sufficiently large values 
of P, individuals become more rational and are able to reason that the 
benefits of taking vaccination can significantly overweigh its cost. In 
this case, many nodes (including some hub nodes) choose to be 
vaccinated. The vaccination coverage reaches a stable level for very 
large values of p. However, for larger values of i5, the vaccination 
coverage decreases with p, as shown in Fig. 2(a3) for i5 = 0.7. Note 
that a large value of S means small cost of vaccination (c(l — S) = 
0.27), so the behaviors are similar to those for c = 0.5. 

For the free-subsidy policy, the vaccination coverage depends 
upon the value of P in a non-monotonic fashion. For small values 
of p, due to the existence of donees (vaccinated without personal 
cost), there is probability that certain individuals imitate these 
donees' strategy even though the vaccination cost is large. As these 
individuals become more rational (larger value of P), they tend to 
switch their vaccination strategy due to the large cost, leading to a 
sharp decline in the vaccination coverage and a rapid increase in the 
epidemic size correspondingly. As the vaccination coverage is 
increased, the risk of infection is reduced, and eventually the cov- 
erage is maintained at a fixed level for large values of p. The insets in 
Fig. 2 show the values of Vp — Vp and of Rp — Rp from numerical 
simulation and theory with reasonable agreement. 

From Figs. 1 and 2, we can observe an interesting phenomenon. 
Specifically, the vaccination coverage for the free-subsidy case is 
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Figure 1 | Quantitative characterization of impacts of the selection strength p on epidemic spreading for partial- and free- subsidy policies for small 
vaccination cost. For homogeneous small-world networks (HSWN) and cost of vaccination c = 0.5, vaccination coverage (V, top panels) and epidemic 
size (R, bottom panels) for the two subsidy policies for c) = 0.1, 0.4 and 0.7 (corresponding to the left, central, and right panels, respectively). Insets in the 
top panels show the difference in the vaccination coverage, Vp — Vp, between the partial- and free-subsidy policies. Insets in the bottom panels display the 
difference in epidemic size, Rp — Rp, between the two subsidy policies. Blue squares in the insets are simulation results, and the red circles are the 
theoretical predictions based on the mean-field method (see analysis in Methods). Other parameters are: network size N = 1000, transmission rate 
A = 0.072, and average degree k = 10. 



larger than that for the partial-subsidy case in most cases, indicating 
that the former policy is more effective at controlling epidemic out- 
breaks. Qualitatively, this phenomenon can understood, as follows. 
In the small fl regime, although a blind imitation behavior exists in 



both cases, the donees in the free-subsidy case never switch their 
strategies, leading to a larger probability for them to allure their 
neighbors into taking vaccination. This can be verified more expli- 
citly by comparing the average fractions of imitated vaccinators 
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Figure 2 | Impacts of the selection strength ^ on partial- and free-subsidy policies for large vaccination cost. For HSWN and c — 0.9, vaccination 
coverage y (top panels) and epidemic size R (bottom panels) for the two subsidy policies for different values of (5 ((5 = 0.1, 0.4 and 0.7, corresponding to 
left, central, and right panels, respectively). Other parameters are the same as for Fig. 1. 
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Figure 3 | The role of donees in epidemic dynamics under free-subsidy policy. Average fractions of imitated vaccinators among neighbors of donees ( ) 
and among non-subsidized vaccinators {Y„i) for c = 0.5 (top panels) and c = 0.9 (bottom panels). Other parameters are the same as for Fig. 1. 



(excluding donees) among donees' and non-subsidized vaccinators' 
neighbors, denoted by and Vnd-, respectively, for the free-subsidy 
case, as shown in Fig. 3. We observe that, for different values of c, the 
value of is no less than that of y„j, suggesting that the donees act 
as the "role models" and are more appealing than other vaccinators 
in attracting their neighbors to take vaccination. To gain further 
insights, we show in Fig. 4 snapshots of the states of individuals in 
a square lattice with periodic boundary conditions for the two pol- 
icies for c = 0.5, fS = 1.0 and S = 0.4, where the vaccination coverage 
(green) in the partial-subsidy case is lower than the free-subsidy case, 
leading to larger epidemic size (yellow) in the former case. As specu- 
lated, for the free-subsidy case [Fig. 4(b)], the freely vaccinated 
donees (red) are randomly distributed in the square lattice and they 
have many surrounding imitated vaccinators, forming vaccination 
clusters, and leading to relatively large vaccination coverage as com- 
pared with the partial-subsidy case [Fig. 4(a)]. 

For intermediate values of /i, individuals possess some degree of 
rationality, so their imitation behaviors are affected by the value of c 
to certain extent. For example, for c = 0.5 (top panels of Fig. 1), as /? is 
increased, while there is a sharp decrease in the vaccination coverage 
for the free-subsidy case, the coverage is still larger than that for the 
partial-subsidy case due to the relatively small value of c. However, 
for c = 0.9 (top panels of Fig. 2), the behaviors are different. 
Especially, when the values of 5 are not too large [Figs. 2(al-a2)], 
the vaccination coverage for the free-subsidy case is lower than that 
for the partial-subsidy case. There are two reasons for this. First, the 
fraction of donees in the population is small, as shown in Fig. 5(a) 
(red and black lines with 5 = OA and 5 = 0.4, respectively). Second, 
large cost discourages other individuals from getting vaccinated. For 
example, from Fig. 5(b), it can be seen that the fractions of the 
imitated vaccinators are about 0.1 for S = Q.l and 0.17 for S = 0.4 



for /? ~ 3.0. For large values of 8, the fraction of donees in the 
population is large, as shown in Fig. 5(a) (blue line). In this case, 
the vaccination coverage for the free-subsidy case is still larger than 
that for the partial-subsidy case even though the fraction of imitated 
vaccinators is very small [blue line in Fig. 5(b)]. These results heur- 
istically explain why, imder the free-subsidy policy, the vaccination 
coverages for d = 0.1 [Fig. 2(al)] and S = 0.4 [Fig. 2(a2)] are smaller 
than that under the partial-subsidy policy but this phenomenon does 
not occur for 8 = 0.7 [Fig. 2(a3)]. 

In the large P regime (e.g., /? s 40 in Figs. 1 and 2), the vaccination 
coverage for the free-subsidy case is again larger than that for the 
partial case, regardless of the values of 8 and c. In this regime, indi- 
viduals are extremely sensitive to the payoff difference. From Eq. (1), 
the probabilities for an individual to switch from a vaccinator to a 




Figure 4 | Typical snapshots of stationary state configuration in a square 
lattice. Panels (a) and (b) are for partial- and free-subsidy policies, 
respectively. Other parameters: network size N = 100 X 100, c = 0.5, 
/? = 1.0 and i5 = 0.4. Blue, green, yellow, and red points represent free- 
riders (not vaccinated and healthy), imitated vaccinators, infected 
individuals, and freely vaccinated donees, respectively. 
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Figure 5 | Fractions of imitated vaccinators and donees under free- 
subsidy policy. For c = 0.9, fractions of donees (a) and imitated 
vaccinators (b) versus /? for different values of (5. Other parameters are the 
same as for Fig. 1 . 

free rider and from being infected to a vaccinator are unity, and the 
other probabilities in the equation are zero insofar as 0 < c < 1. As a 
result, the vaccination coverages for the two subsidy cases are not 
dependent upon the value of c any more. Consequently, the coverage 
for the partial-subsidy case no longer changes with the value of b 
either, because of the condition 0 < c(l — (5) < 1. (Numerical 
support for this can be found in Figs. 1 and 2, where the vaccination 
coverage for the partial-subsidy case is about 0.37.) Under the free- 
subsidy policy, if there are no donees in the population (i.e., p = 0), 
the vaccination coverage should be approximately the same as that 
for the partial-subsidy case. Furthermore, for the free-subsidy case, 
since some donees never switch their strategies and thus are capable 
of attracting some of their neighbors to take vaccination, a higher 
level of vaccination coverage would arise. This reasoning can also be 
verified theoretically (see Methods). 

For even larger values of [i, the equilibrium solution x in Eqs. (18) 
and X in Eq. (20) can be obtained by using Eqs. (12) [or Eq. (17)] and 
(19). In particular, Eq. (18) shows that the vaccination coverage for 
the partial-subsidy case is not dependent upon the value of c or S. 
That is, when individuals become completely rational, the partial- 
subsidy policy has little effect on the epidemic dynamics. 
Nevertheless, for the free-subsidy case, Eq. (20) stipulates that the 
vaccination coverage is larger and in fact increases with the fraction 
of donees, p, where p = 8 • Vp, with Vp being the vaccination cov- 
erage for the partial-subsidy case. This gives theoretical support for 
the observation that the donees under free-subsidy policy can 
enhance the level of vaccination when individuals are completely 
rational. Figure 6 shows the simulation results for fi = 50, together 
with the theoretical prediction from Eqs. (18) and (20). It can be seen 
that the vaccination coverage for the partial-subsidy case is nearly 
constant but the coverage for the free-subsidy case increases with the 
value of (5. 

Figures 7 and 8 show the results from mean-field theory 
(Methods) for a systematic comparison of the effects on epidemic 
dynamics of the two distinct subsidy policies. In particular. Fig. 7 
compares the effects of the two subsidy policies with respect to var- 
iations in the parameter plane {fS-8) for c = 0.5 and c = 0.9. From 
Fig. 7(al), one can observe that, for c = 0.5, the vaccination coverage 
for the partial-subsidy policy is no larger than that for the free-sub- 
sidy case, rendering inefficient control of outbreak of epidemics 
using the partial-subsidy policy [Fig. 7(bl)]. For c = 0.9, the vaccina- 
tion coverage for the partial-subsidy case can exceed that for the free- 
subsidy case but only for the small S regime, as shown in Figs. 7(a2) 
and 7(b2). In Fig. 8, the comparison is carried out in the p-c para- 
meter plane for (5 = 0.1 and S = 0.7, where it can be seen that the 
vaccination coverage for the partial-subsidy case can exceed that for 
the free-subsidy case only for c close to unity. However, the difference 
diminishes when S becomes large, as shown in Fig. 8(b). Taken 
together, these results suggest that the free-subsidy policy can in 
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Figure 6 | Stable vaccination coverages under free- and partial-subsidy 
policy for very large value of p. For /? = 50, numerically obtained 
difference in the vaccination coverages between the two cases, Vp — Vp, 
and the corresponding theoretical prediction from Eqs. (16)-(18). Other 
parameters are the same as for Fig. 1 . 

general control epidemic outbreak more effectively than the par- 
tial-subsidy policy, except for large values of c and intermediate 
values of p. 

Real complex networked systems often possess certain degree of 
skewness in their degree distributions, typically represented by some 
scale-free topology. We thus implement our model of human-beha- 
vior based epidemic dynamics on the Barabasi- Albert (BA) scale- 
free" networks to assess the generality of the phenomena uncovered 
for regular and homogeneous networks. Figures 9 and 10 illustrate 
the effects of selection strength [i under the two subsidy policies for 
c = 0.5 and c = 0.9 (c = 0.2 is given in Fig. S2.), respectively. (Results 
from the configuration network'^ are also given in Figs. S3-S5.) We 
observe essentially similar results to those for HSWN networks 
(Fig. 1 and Fig. 2). Theoretical results under the two subsidy policies 
based on Eqs. (3)-(9) are given in the insets of Figs. 9 and 10. We 
observe a reasonable agreement between theory and numerics. 

Discussion 

The spontaneous behavioral responses of individuals to epidemic 
spreading have a significant impact on its dynamics. In the past 
decade, there was a great deal of interest in the development of 
mathematical models attempting to close the epidemic spreading 
behavioral changes epidemic-spreading feed-back loop^'"'". 
In the real world, various subsidy policies are often devised to 
encourage the individuals to actively participate in disease-preven- 
tion programs. However, governmental incentive programs can alter 
the behavioral responses of individuals, reducing the effectiveness of 
these programs and even leading to detrimental consequences. 
Studies of epidemic models should then take into account the inter- 
play among spreading dynamics, individual behaviors, and public 
policies to yield insights into the complex dynamical process with 
useful and predictive results. 

We propose a vaccination-decision model based on evolutionary 
game theory to assess the effectiveness of two realistic subsidy pol- 
icies, the free- and partial-subsidy policies, in a quantitative manner. 
In particular, using a physically reasonable parameter, the selection 
strength, we uncover the complex and distinct roles played by the two 
policies in controlling epidemic spreading. When the partial-subsidy 
policy is imposed, for small cost of vaccination, individuals' being 
more rational can lead to smaller vaccination coverage. For high 
vaccination cost, as the individuals become more rational the vac- 
cination coverage tends to increase. Under the free-subsidy policy, at 
small vaccination cost the eventual extent of vaccination exhibits 
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Figure 7 | Comparison between the effects of free- and partial-subsidy policies in the parameter plane (p, rf). Left and right panels: differences in the 
vaccination coverage and epidemic size, Vp — V^and Rp — Rp, respectively, in the pS plane for c = 0.5 (al-bl) and c = 0.9 (a2-b2). Other parameters are 
the same as those in Fig. 1 . 
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Figure 8 | Comparison between the effects of free- and partial-subsidy policies in the parameter plane {p, c). Left and right panels: differences in the 
vaccination coverage and epidemic size, Vp — V^and Rp — Rp, respectively, in the /?-cplanefor (5 = 0.1 (al-bl) and (5 = 0.7 (a2-b2). Other parameters are 
the same as for Fig. 1. 
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Figure 9 | Impacts of tlie selection strength ^ on epidemic dynamics for scale-free networks for small vaccination cost. Vaccination coverage V (top 
panels) and epidemic size R (bottom panels) under the two subsidy policies for scale-free networks for c = 0.5 and d = 0.1, 0.4 and 0.7 (corresponding to 
the left, central, and right panels, respectively) . Insets in top and bottom panels show the difference in the vaccination coverage, Vp — , and the 
difference in the final epidemic size, Rp— Rp, respectively, between the two subsidy policies, where the blue squares and red circles are simulation results 
and mean-field based theoretical prediction (Methods), respectively. Other parameters are: N = 1000, / = 0.18, and average degree (k) = 6. 



similar behavior with respect to variations in the individual ration- 
ality to that under the partial-subsidy policy. However, for high 
vaccination cost, a complex chain of relationships emerges. 
Specifically, for a population with relatively less rational individuals, 
high vaccination coverage can arise. Having a more rational popu- 
lation of individuals can result in a sharp decline in the vaccination 



coverage. For a population of extremely rational individuals, vaccina- 
tion can be substantially favored again. This counterintuitive phe- 
nomenon, namely, under the free-subsidy polity relatively high 
vaccination coverages can occur in a population of either very irra- 
tional or very rational individuals, can be explained by focusing on a 
natural type of behavioral responses observed commonly in human 
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Figure 10 | Impacts of the selection strength ^ on epidemic dynamics under the partial- and free-subsidy policies for scale-free networks for large 
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societies: imitation. For irrational individuals, blind imitation beha- 
viors promote the role models of the donees under the free-subsidy 
policy, leading to a higher vaccination coverage. For very rational 
individuals, the existence of the donees leads to a higher stable equi- 
librium in the vaccination coverage under the free-subsidy than 
under the partial- subsidy policy. All these results have been estab- 
lished by theoretical analysis and extensive computation for both 
homogeneous and heterogeneous complex networks, with good 
agreement between the respective results. Taken together, a free- 
subsidy strategy to control disease can be more effective than a par- 
tial-subsidy strategy, but the opposite holds if the level of rationality 
in the society is neither too low nor too high. 

Our efforts may help push one step forward the field of quantitat- 
ive analysis and mathematical modeling of epidemic spreading and 
disease control, an problem of paramount importance and broad 
interest in the modern time where rapid epidemic spreading on a 
global scale represents one of the most dangerous threats to the 
human race. Our findings reveal that counterintuitive phenomena 
can arise due to the complex interplay among epidemic dynamics, 
human behaviors, and governmental policies. Especially, the making 
and execution of any disease- control policy should be exercised with 
extreme caution: its success depends not only on a number of rel- 
evant factors but, more importantly, on their dynamical interplay. 
Our work provides considerable insights into the possible conse- 
quences of such interplay, calling for further efforts in this interdis- 
ciplinary field. 

Methods 

Model. Taking into account periodic outbreaks of flu-like diseases and the limited 
effectiveness of vaccines, we study models with two dynamical ingredients: seasonal 
updating and pre-emptive vaccination, in which individuals decide whether or not to 
get vaccinated before each epidemic season. Based on a previously studied modeP *^, 
we integrate the two dynamical ingredients by generating an iterative system of a two- 
stage process. Firstly, individuals make decisions during a yearly vaccination 
campaign - the vaccination stage. A vaccinated individual pays a cost Cy that depends 
on issues such as the time needed to get vaccinated (e.g., via public health services) 
and potential side effects*^"*^. For simplicity, we assume that the vaccine provides full 
protection against the infection in the next epidemic season. The second stage is the 
epidemic season, during which each susceptible and non-vaccinated individual has 
certain probability to become infected sometime during the season. In this stage, the 
epidemic strain infects an initial number Iq of individuals, and the disease spreads 
according to susceptible-infected-resistant/vaccinated (SIR) dynamics, with daily 
transmission rate A and recovery rate //. The epidemic continues until there are no 
more newly infected individuals. An individual who gets infected during the epidemic 
season pays a cost Q. Without loss of generality, we rescale the payoff by defining the 
relative cost of vaccination as c — CylCf. Generally, the vaccination cost should be 
smaller than the cost of infection so that 0 < c < 1; otherwise, there is no incentive to 
be vaccinated. Those individuals who are neither vaccinated nor infected pay no cost, 
so they are the free riders. 

Once the epidemic ends, individuals update their vaccination decisions for the next 
season by imitating the strategy of neighbors who have gained higher payoff. A 
commonly used updating strategy in evolutionary- game theory is one given by the 
Fermi rule. By this rule, an individual, say i, updates his/her vaccination strategy by 
randomly choosing one of the immediate neighbors, say j, compares their costs (or 
payoff), and adopts the strategy oij with the following probability as determined by 
the payoff difference^^: 



l + exp[-/?(P,--P,)] ■ 



(1) 



where s, — 1 or 0 denotes the vaccination choice for individual i: either vaccinated or 
not, and Pj is the current payoff of individual i at season t. Without any subsidy, 
according to the costs of vaccination and infection, we have 



— c, vaccination; 

— 1, infected; 
0, free-rider. 



(2) 



where ^ is a parameter characterizing the rationality of the individuals. In particular, 
irrational and very rational individuals correspond to and ^;^>1, respectively. 

We implement the above epidemiological/game-theoretic model on HSWNs, 
constructed by random link rewiring from a regular network such that the degree of 
each node remains unchanged. In particular, a pair of edges, A — B and C — D, are 
randomly selected. They are then rewired to generate new link pairs A — D and B — C. 
Multiple edges connecting the same pair of nodes are prohibited. The process 



continues for PE steps, where E is the number of edges in the network and P char- 
acterizes the randomness of the network. For P > 1 the network topology is essentially 
random^**. (To be concrete, we set P — 5.O.). We also use the BA scale-free networks^\ 
In our simulations, the initial state is composed of equal fractions of vaccinated and 
non- vaccinated individuals, who are randomly distributed throughout the entire 
population. The number of initially infected individuals is Iq — 5. We calibrate the 
value of disease transmission probability to ensure that the infection risk in an non- 
vaccinated population is equal across all possible contact configuration. The trans- 
mission rate is chosen to be A — 0.072 per day per person for HSWN with average 
degree 10, and a — 0.18 per day per person for BA scale-free network with average 
degree 6. The value of // is chosen to be — 0.25 per day for all contact pairs to ensure 
that the fmal infection size is about 0.9^'^^. We carry our 100 independent statistical 
realizations and, for each realization, the simulation is run for 3000 iterations, where 
the last 1000 iterations are used to calculate the statistical quantities about the epidemic. 

Steps are also taken to ensure fair comparison between the two types of subsidy 
policies. In particular, for partial-subsidy policy, each vaccinated individual has his/ 
her vaccination cost reduced by the relative amount 3, so the actual cost paid by the 

individual is i:(l — 3). The total amount of subsidy is c3 ^v{t), where Tis the 

total number of iterations and Ny (t) is the total number of vaccinated individuals at 
season t. Under the free-subsidy policy, the total amount of subsidy is equally allo- 
cated to each vaccination stage, so the number of freely vaccinated donees at each 

stage is c3 ^'^^^^vit) /{Tc)^3 ^'^^ ^ Nv{t) /T. Since, in this case, the number of 

donees at each vaccination stage is unchanged, we do not choose donees randomly at 
the following vaccination stages. That is, if one individual is selected as a donee for the 
initial season t — 1, then the individual will be a donee for subsequent seasons. 

Analysis of epidemic dynamics in networks. For a complex network with degree 
distribution P{k), we let S^it), hit), Rkit) be the densities of susceptible, infected and 
recovery nodes of degree k at time t, respectively, where S^(f) + Ii^{t) + Rk{t) = 1. 
Under the mean-field approximation, the time evolutions of the densities are given by 
the following set of coupled differential equations^": 



dSiJj) 

dt 

dh{t) 

dt 
dRkit) 

dt 



= ).kS^(t)&(t)-nh{t), 



(3) 



= tih{t). 



where the factor Q{t) — k'P{k'\k)Ik'{t) represents the probabihty that any given 
hnk points to an infected node. In absence of any degree correlations, Q{t) is given 
by=' 

0(t)=^m'c)4(o/('^)- 



Let xi^ be the density of vaccinated nodes with degree k before the epidemic season. 
We have 



dh{t) 

dt 
dRkjt) 
dt 



--Xk{l-xt-h{t) -Rt(t))&{t)-tilt{t), 



(4) 



After each epidemic season, one can obtain the degree-dependent fraction of epi- 
demic size Rk{^). 

Without any subsidy policy, the payoff of each individual after each SIR epidemic 
season is given by Eq. (2). Specifically, for a vaccinated individual (V), the payoff is Py 
— —c. For an non -vaccinated and healthy individual (H), we have P/j — 0. For an 
non-vaccinated and infected individual (7), the payoff is P/ — — 1. Individuals are 
allowed to switch their vaccination strategies with probability given by Eq. (1). 

Under the partial- subsidy policy, a vaccinated individual's payoff is modified as 
Psv — (1 ~ where SV stands for subsidized vaccinators. In this case, whenever an 
SV changes to H or / state, the variable Xj^ is decreased according to 



(5) 



where the first and second terms on the right-hand side of Eq. (5) denote the tran- 
sitions from vaccinators to infected individuals and from vaccinators to free-riders, 
respectively. In the opposite case where an individual in H or / state is changed to 
being an Sy individual, is increased according to 



^/cP(/c)XfcWj 



(6) 
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Combining Eqs. (5) and (6), we can write the time evolution of tlie vaccination 
coverage as 



(7) 



Here, to avoid confusion, we use t to specify the epidemic season and t as the time 
variable for epidemic seasons. The equilibrium vaccination coverage and the fmal 
epidemic size can be calculated by numerically solving Eqs. (4) and (7) iteratively. 

On average, the total amount of subsidy at each vaccination stage is NxcS for the 
partial- subsidy case with x — P{k)xk, so the fraction of freely vaccinated donees 
under the free-subsidy policy isp = NxcS/{Nc) — xS. Similar to Eq. (5), the reduction 
in xj^ , with freely vaccinated donees included, can be written as: 

x^^'^^i^J^^Pik^Rki'^^^^ (8) 

where xi^ — p is the fraction of the individuals with degree k that can change their 
vaccination strategy. Meanwhile, it is not necessary to replace p by since the donees 
are randomly selected. The gain of for the free-subsidy case is given by 



{1-Rk{x,)-Xk) 



J2kP(k}xtWH^V. 



(9) 



For HSWNs under the partial- sub sidy policy, Rki^) and x^ in Eqs. (5) and (6) can be 
replaced by andx, respectively. By using Eq. (13) below [R{^) — (1 — x)f{x)], we 
can simplify Eqs. (5) and (6) as 



x-^x{l-x) \f{x) Wsv^i + ( 1 -/ W ) Wsv^h] , 



and 



x+^x{l-x) lf{x) Wj^sv + ( 1 -/ W ) Wh^sv] . 
After some straightforward algebra, we have 



dx{x) 
dz 



'{ 



= x{l-x)i [1 -/(;>:)] tanh 



+/(z)tanh||(Psv,-P/) 



(10) 



(11) 



(12) 



= x( 1 - x) I [1 -fix)] tanh ( - c( 1 - ^))] +f(x) tanh ( - <:( 1 - a) + 1)] } , 

where /(a:) is the probability that a susceptible individual finally gets infected in a 
population with vaccine coverage x, which is given by 



R(oo) 
1-x ' 



The final epidemic size, R{^), can then be determined by 

dl{t) 



dt 

dR(t) 
dt 



= Ak{l-x-I{t)-R{t))I{t)-fiI{t), 



(13) 



(14) 



with k being the degree of each node. Using a similar method of analysis for well- 
mixed network^-^, we obtain the following self-consistent equation from Eq. (14): 



R(oc.)-(l-jc) l-e-«'^^) , 



15 



where i^o — kXj ^iis the basic reproduction ratio of the epidemic. Combining Eqs. (13) 
and (15), we have 



x=\ + 



l"(l-/(^)) 
«o/(^) 



Setting the right-hand side of Eq. (12) to be 0, we get 



/(^) 



tanh 


(^<i 


-)] 




tanh 


^^^(1-^)) 


- tanh 


'^-{[c(\-S)-\)) 



(16) 



(17) 



Substituting Eq. (17) into Eq. (16), we can obtain the steady vaccination coverage. 

For large values of ^, most values of c lie in the interval 1/^ < (1 — b)c < 1 — 1/^ 
since 1/^ ^ 0. In this case, we have 



/(x) = l/2|tanh(^c(l-^) 



- 1 and tanh 



(c(l-^)-l) 



1 . Then, 



by substituting/(x) — 1/2 into Eq. (16) we have the following equilibrium solution 
under the partial-subsidy policy; 



21n2 



(18) 



For HSWNs under the free-subsidy policy, by simplifying Eqs. (8) and (9) and writing 
them as Eq. (12), we have 



dx(x) 
dx 



= x(\-x)\(\-f{x))Wn^y +f{x) W,^v] 

-(x-p){l-x)l{l-f (x))Wy^H +f {X) Wy^,] . 



(19) 



We can obtain the formula of/(jJc) by letting the right-hand side of Eq. (19) be 0, and 
then substituting it into Eq. (16). 

The steady-state vaccination coverage for the free-subsidy case can be numerically 
solved too. Especially, for large values of ^, and c satisfying 1/^ < c < 1 — 1//5, the 
switch probabilities in Eq. (19) can be approximated as: Wh^v — 0, Wj^y ~ 1> 
Wv-*H ~ 1 Slid Wy^i — 0, so we have/(x) — {x — p)l{2x — p). Substituting it into Eq. 
(16), we obtain the following self-consistent equation for the equilibrium value x for 
the free-subsidy case as 



x^l + - 



{ 2x-p)ln^^ 
Rq(x-p) 



(20) 
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